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S ' Abstract 

\ We address the question of cosmological perturbations in the context of 

brane cosmology, where our Universe is a three-brane where matter is con- 
psj ■ fined, whereas gravity lives in a higher dimensional spacetime. The equations 

governing the bulk perturbations are computed in the case of a general warped 
^ , universe. The results are then specialized to the case of a five-dimensional 

\ spacetime, scenario which has recently attracted a lot of attention. In this 

^ I context, we decompose the perturbations into 'scalar', 'vector' and 'tensor' 

VO ' modes, which are familiar in the standard theory of cosmological perturba- 

tions. The junction conditions, which relate the metric perturbations to the 



o 
o 



, matter perturbations in the brane, are then computed. 

> , 

^ I. INTRODUCTION 

D ; 

Recently has emerged the fascinating idea that extra dimensions (beyond our familiar 
time and three spatial dimensions) could be large today or at energy scales much below what 
^ \ was thought before The reason why these extra dimensions could remain hidden is that 
d ' matter fields would be confined to a three-dimenional brane (our visible 'space') whereas 
gravitational fields would live everywhere and thus also in the extra- dimensions. 

Even non-compact extra-dimensions can be envisaged 0, like in the five-dimensional 
model of Randall and Sundrum 0, where the bulk is endowed with a negative cosmological 
constant. Remarkably, in this model, the gravity felt by observers on the brane will behave, 
at least at first approximation, like usual four- dimensional gravity, because of the existence 
of zero-mode gravitons effectively confined in the brane. 

When one considers a brane universe in the context of cosmology, it appears that one 
does not recover the standard Friedmann equations easily IQ. The reason is that the matter 
content of the brane enters quadratically in the equations governing the dynamics of the 
brane geometry. However, when one generalizes the idea of Randall and Sundrum to cos- 
mology, i.e. when one assumes the existence of a negative cosmological constant in the bulk, 
adjusted so as to compensate the tension of the brane, one ends up with a cosmological 
evolution that is conventional at sufficiently low energies ( and [^]). 

The question now is whether this conventional behaviour will still be recovered in brane 
cosmology, when one relaxes the assumption of homogeneity and isotropy in the brane, in 



1 



other words when one considers inhomogeneous brane cosmology. Up to now, almost all 
works (see e.g. 0) on the cosmological aspects of the brane scenario have dealt with a ho- 
mogeneous and isotropic brane universe. However, it is well known that a lot of information 
is contained in the cosmological perturbations of our Universe, and it is therefore of the up- 
permost importance to analyse the behaviour of cosmological perturbations in the context 
of brane cosmology. It would be crucial to test whether brane cosmology is compatible with 
current observations of cosmological perturbations, either the anisotropies of the Cosmic 
Microwave Background or the large scale structure data. It would also be interesting to 
check if the usual mechanism of generation of cosmological perturbations via amplification 
of quantum fluctuations during an inflationary phase would still be valid in brane cosmology 
(see i). 

Another motivation to study the perturbations in the cosmological context is the question 
of the stabilisation of the radion and its implication on the cosmological evolution in the 
brane P|. As shown by |]TU[ in a non-cosmological context, a perturbative approach may be 



necessary for a full understanding of the radion. 

The purpose of the present paper is to develop a formalism that can describe the evo- 
lution of the cosmological perturbations in the context of brane cosmology. In order to 
do so, we will proceed in several steps. First, one must compute the perturbations of the 
bulk Einstein's equations, which relate the perturbations of the metric to the perturbations 
of the bulk matter, if there is any. For this first step, we have computed the perturbed 
Einstein's equations in the more general case of a warped spacetime with any number of 
dimensions. A similar calculation has been done recently in the case of maximally 
symmetric spacetimes. 

The results valid for any dimension can be applied to the model we wish to focus on: 
a five- dimensional spacetime. It is then useful to distinguish several types of cosmological 
modes, following the familiar decomposition of the standard theory of cosmological pertur- 
bations (see e.g. or [|T^). 



At this stage, however, the perturbations of the matter in the brane have not yet been 
taken into account. They play a role only in the junctions conditions: the jump in the 
derivatives (with respect to the fifth dimension) of the bulk metric perturbations is indeed 
governed by the matter perturbations in the brane. This is in this very unusual way that 
matter perturbations in our apparent Universe, i.e. the three-brane, are connected with 
the geometry perturbations in our Universe, i.e. simply the particular value on the brane 
of the bulk metric perturbations. It is nevertheless interesting to notice that the present 
problem has some analogy with the question of the interaction between a domain wall and 
gravitational waves in a four-dimensional spacetime 

The plan of the paper will be the following. In Section 2, we obtain the perturbations 
of Einstein's equation in the bulk, in the case of a general warped spacetime. In Section 3, 
we consider the five- dimensional model and recall what is known about the homogeneous 
and isotropic solutions in an empty bulk or with a negative cosmological constant. Section 
4 is devoted to the bulk perturbations of the five-dimensional spacetime, whereas section 5 
deals with the junction conditions that take into account the matter in the brane. Finally, 
we conclude in the last section. 
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II. BULK PERTURBATIONS IN GENERAL WARPED SPACETIME 



Although we will be ultimately interested in the perturbations of a five- dimensional 
spacetime, it is instructive to begin with a more general situation and to compute the 
perturbations of the Ricci tensor for a D-dimensional spacetime that can be considered as a 
warped product of a p-dimensional spacetime with a d-dimensional space. The metric has 
thus the particular form 

gAsdx'^dx^ = ^abdx"'dx'^ + a'^{x'^}'yijdx''dx^ , (1) 

where A, B, ... denote global spacetime indices, a,b, ... indices of the p-dimensional spacetime 
with metric 'jab , and ... indices of the d-dimensional space with metric 'jij. 

Let us denote the global covariant derivative associated with the metric (Jab, whereas 
Va will stand for the covariant derivative associated with the metric jat and Vj for the 
covariant derivative associated with the metric 7jj. The nonvanishing mixed Christoffel 
symbols (those with indices of both types a and i) are 

^a, = ^lh ^ = -{acfaH,. (2) 

All tensors components will be decomposed in several sets, depending on the number of 
p-indices and of d-indices. And one can then decompose the action of the global spacetime 
covariant derivative Da into the action of the covariant derivative Va and that of the covari- 
ant derivative Vi. For illustration, on a vector = (y^^V^), the action of the covariant 
derivative Da gives 

DaV = VaV\ 

DiV = ViV - iad''a)-fijV\ 

D,v^ = + ^7;.r^ (3) 

This type of formulas can be generalized to any type of tensor by using the mixed Christoffel 
symbols given above. 

Let us now consider a linear perturbation of the spacetime metric, so that the total 
metric reads 

ds^ = {jjAB + hAB)dx^dx^ = gABdx^dx^ . (4) 

Our purpose will now be to compute the linear perturbation of the Ricci tensor, 6Rab- 
Quite generally, the expression of 6Rab in terms of the metric perturbations is found to be 
of the following form (see e.g. |T5[ ) 

5Rab = —DADsh - ^D^DchAB + D^'D^BhA)c, (5) 
where h denotes the trace of the perturbation Hab, namely 

h = g^'^'hAB. (6) 
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After some tedious but straightforward calculations, it is possible to express the various 
components of the linearized Ricci tensor more specifically in terms of the components hab, 
hai and hij, and of the covariant derivatives Va and Vj. One finds 

+ 2 iyahbc + ^bhac - ^chab) - a~^V i^aO^ b)h, 

1 2^k^ , lV„a^,, ^aa^h, dV''a~, / ci\ V^a ~ , 
+-a-2v'=Vi/i„fc + -^Vi{h -h)- -^^''hbi + Wahbi + 1 - 7T ^bhai 



7^ ^iKb Kk -{d- 1) — hbi } 

2 1a a"* a 



1 1 ~ ~ 1 1 1 

-^V"V„/i,,- + U/'^ViKi + ^V"V,-/i„i + ^lijV^hak + (aV^a) li^V'^h^ 
d\V''a~ , (d A V^a^ , (d \ V^a^ , 

+ {d - l)(V"aV'a)7,,/i,fe + a(V"V'a)7,,/i„b - 2^^^/i,„ 



with ^ = a ^7*-'7i, 



Gauge transformations. 

It is important to notice that perturbations which differ quantitatively can in fact describe 
the same geometry simply because they correspond to different systems of coordinates. 
Gauge transformations, corresponding to infinitesimal changes of coordinates 

x^^x^ + C^, (10) 

induce the following transformations for the metric perturbations, 

hAB^hAB + DA^B + DBU- (H) 

This general expression gives for the three types of metric components or, decomposing into 
the two subsystems of coordinates, 

hob Kb + Ibc^aC + lac^bi", 

hij hij + a' (ijk^ii'' + lik'^ji^) + 2a(9„a)7i,r, 

hia hia + 7a&ViC'' + a^-fijV a^' ■ (12) 



III. FIVE-DIMENSIONAL BACKGROUND SPACETIME 



In this section, we recall briefly the cosmological solutions that have been found in the 
case of five-dimensional spacetimes, with a metric of the form 



ds^ = —n^ir, y)dT'^ + a^(r, y)jijdx^dx^ + b'^^r, y)dy'^, (13) 

where the spatial three-surfaces, defined by r and y constant, are homogeneous and isotropic 
and is a maximally symmetric 3-dimensional metric {k = —1,0,1 will parametrize the 
spatial curvature). 

The five- dimensional Einstein equations take the usual form 

Gab = Rab — -^Rqab = i^^Tab-, (14) 

where Tab is the five- dimensional energy-momentum tensor. 

With the above metric, the non-vanishing components of the Einstein tensor Gab are 
found to be 

G =3/^ + (^^ + ^ +A;— 1 (15) 

°° [aya hj IP' \a a \ a h j j a^J' 

a" ia' fa' n'\ h' ( n' a'\ a" n"\ 

(r ya\a n J b \n a J a n j 

a? { a / d h\ d b / d h\ b] , , ^ 

[a (" + ^n) - n ("^a + n) - l] - '^^^ 

G,,^i{'-^'^+'^\-'^]. (17) 

\n a a b a J 

^ o n'\ b"^ fd fd n\ d\ h"^] 

^55 = 3 - - + - -— +- -A:- , 18 

[aya n J n^VaVa n/ a/ a"^J 

where a prime stands for a derivative with respect to y, and a dot for a derivative with 
respect to r. The stress-energy-momentum tensor can be decomposed into two parts, 

where ^Luik the energy momentum tensor of the bulk matter and sLrano corresponds 
to the matter content in the brane (y = 0). Since we consider in this section only strictly 
homogeneous and isotropic geometries inside the brane, the latter will be necessary of the 
form 

T^bL... = ^diag(-p,p,p,p,0), (20) 

where the energy density p and pressure p are functions only of time. 

In the case of a general stress-energy tensor for the matter on the brane, and assuming 
a time-independent metric along the fifth dimension, explicit solutions for the whole metric 
have been found in two simple cases: the case where the bulk is empty @] and the case 
where the bulk contains a cosmological constant ||^, i.e. 

slbuik = diag {-pB, -pB, -pB, -pB, -Pb) , (21) 

with Pb = const. 
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In the case of a negative bulk cosmological constant, the most reahstic, the exphcit 
expressions for the metric components (with 6=1 and n{t, y = 0) = 1) are given by ^ 



a{t,y) 




3C 



6p 



B 



cosh(/z?/) 



(22) 



with p 



-"^Pbi and 



n{t,y) 



a{t,y) 



(23) 



where Oq is the scale factor in the brane (i.e. in y = 0). The two functions ao{t) and p{t), 
which appear in the above solution, are obtained by solving 



C k 

an 6 36^ ttn ttn ' 



(24) 



where C is a constant (of integration), and 



p + 3—{p + p) 
ao 



0. 



(25) 



The last equation corresponds to the ordinary energy conservation law. But (24), analogous 
to the first Friedman equation, does not give the usual cosmological evolution because of the 
quadratic term p^. It is however possible to recover the standard cosmological evolution, at 
least at sufficiently late times, if one decomposes p into a tension a and an energy density 
of ordinary matter living in the brane and if one assumes a negative cosmological constant 
in the bulk that will compensate the o"^ term (see [H). 



IV. PERTURBATIONS OF A FIVE-DIMENSIONAL SPACETIME 

After having described our reference spacetime in the previous section, we will now study 
linear perturbations about it. For this task, the calculations of Section 2, specialized to the 
case of five-dimensional spacetime, are going to be very useful. As it is clear from the 
previous section, the metric %b will be defined by 

^abdx'^dx^ = -n{t, yfdt^ + h{t, yfdy^. (26) 

The second metric jij is simply the metric covering the three ordinary spatial dimensions. 
In a cosmological context, there are three choices for jij depending on the spatial curvature 
of spacetime. For simplicity, it will be assumed that our Universe is spatially fiat, which 
means that 

lij = ^ij- (27) 
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From now on, we also choose to work in a Gaussian normal (GN) system of coordinates 
adapted to the three-brane, in which the brane is localized in y = and the metric has the 
form 

ds'^ = g^ydx^'dx'' + dy'^, (28) 

where the greek indices refer to ordinary spacetime coordinates, i.e. the time coordinate 

r and the three ordinary spatial dimensions x*. For the background, this choice of gauge 
corresponds simply to set 

h{t,y) = l. (29) 

In the case of the perturbed spacetime, described by the metric (^, the choice of a GN 
system of coordinates will impose 

/i55 = h^^. = 0. (30) 

Taking this into account, the linearized metric (^) can now be written, quite generally, 
in the form 

ds"^ = -n^(l + 2A)dt^ + 2Bidx'dt + c? (bij + hi^ dx'dx^ + dy^. (31) 



Following Bardeen |jT6[ (see also [jT2[)5 the linearized quantities specified above can be de- 
composed further, into so-called scalar, vector and tensor quantities, according to the ex- 
pressions, 

Bi = d,B + Bi, (32) 

where Bi satisfies diB^ = 0, and 

hi, = 2C6ij + 2didjE + 2d(^^Ej) + Eij , (33) 

where Eij is transverse traceless, i.e. diE^^ = 0, E^ = 0. The quantities A, C, B and E are 
usually refered to as 'scalar' perturbations, Bi and Ei as 'vector' perturbations, and Eij as 
'tensor' perturbations. 



A. Gauge transformation 

As explained earlier, the perturbations defined above can be quantitatively different but 
describe the same geometry if they are written in different coordinate systems. In order 
to distinguish gauge effects and physical degrees of freedom, it is useful to write down the 
effect of a coordinate change (or gauge transformation) on all the metric perturbations 
defined above. These transformations follow directly from the general expressions given in 
(0). 

We will parametrize the infinitesimal coordinate transformation by the vector = 
{C,^,Cj^^)- Let us first consider the components /iss, /iqs and hi^, which vanish in a GN 
coordinate system. They transform according to the laws. 
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^55 ^55 + , 

hi5 ^ + di^ + aHi^e'. (34) 

In order to bring any coordinate system into a GN coordinate system, it is clear from the 
above relations that once one has used to adjust the position of the brane at ?/ = 0, 
appropriate choices of and of ^* will be required in order to make /iqs and hi^ vanish. 
Note however that the GN gauge adapted to the brane is not completely fixed: there is 
some residual gauge freedom associated with parameters and that depend only on the 
four ordinary spacetime coordinates. This residual gauge freedom can be interpreted as 
possible redefinitions of the coodinates inside the brane worldsheet. 

Let us now consider the transformations for the other components of the metric pertur- 
bations. To do that, it is convenient to decompose the spatial vector into 

e = a^e+e (35) 

(where (9* = 6^^dj), such that is transverse, i.e. di^^ = 0. With this decomposition, one can 
see that the three scalar parameters ^ and will induce transformations in the subset of 
scalar perturbations, whereas will act in the 'vector' subspace. The tensor perturbations 
Eij will be left untouched by the gauge transformations. Specializing the expressions (|12D 
to our particular case here, one finds the following transformations: 



1. Scalar gauge transformations 
n n 

a a 

E~^E + i, (36) 



2. Vector gauge transformations 
Bi^B, + a%, 

E,-^E, + li. (37) 

In the above transformations, we have considered the most general gauge transformation. 
But if one considers only gauge transformations inside the subset of GN coordinate systems, 
which will be the case in the following, then disappears from (|36|) and the parameters 
^ and ^ cannot depend on the fifth coordinate. 

Finally let us remark that one can use the remaining gauge freedom within the GN 
subset to impose some additional gauge conditions. For example, for scalar quantities, a 
choice which is familiar in the standard theory of cosmological perturbations is to impose 

5 = 0, E = 0, (38) 

which corresponds to the so-called longitudinal (or Newtonian) gauge. However, these con- 
ditions can be imposed only on one hypersurface, y = say, but not everywhere in the bulk, 
because B and E a priori depend on the coordinate y. 
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B. Perturbed Einstein equations in the bulk 



Let us now compute the equations governing the perturbations in the bulk, ignoring for 
the moment the presence of the brane. The five dimensional Einstein's equations ([T^), in 
the bulk, can be rewritten 

Rab = (^Tab - ^Tqab^ ■ (39) 
Therefore, the perturbed Einstein equations in the bulk have the form 

5Rab = {sTab - \ThAB - \9AB8f^ ■ (40) 

Since the background solutions have been found explicity, as recalled in Section 3, the 
cases of an empty bulk or of a cosmo logical constant are of particular interest, but of course, 
the present formalism applies to any bulk energy-momentum tensor. In the case of an empty 
bulk, the perturbed Einstein's equations are simply 

5Rab = 0, (41) 

whereas in the case of a bulk with a cosmological constant A = k^pb, the perturbed Einstein's 
equations read 

SRab = l^hAB- (42) 

The remaining task consists in computing explicitly the components of the perturbed 
Ricci tensor. Using the expressions (0-|^), one obtains the following expressions, which can 
conveniently be separated in scalar, vector and tensor parts. 



1. Scalar components 



SRt 



n 









\ a 


n J 



n 



n 



a'n' 



an 



-3C + 3-A+ 3 



n 



n 



6- C + — AA 
a J 



n 



n 



a 



2-]AE + nn'AE' + a'^AB 



n 



n 



n 



AB, 



n 



a 



n'^ \a 



—a 



' ah a 

4-^-2 — - + 2 — , 



A + 2 



a'n' 



a a n 

+ -A'+ 6- + - 
a \ a n 



an-' 



c 



an^ 
a 



a" 

+ 2— + 
an 



ha 
na 




/ an 


'! 


a" 

— + 
a 


ah 




a 


an^ 


a^n^ 


an^ 



(43) 
(44) 



C 



+C" + — - 6 



n 



a 



an^ 



an^ 
1 



-A 



n^ 



n^d^ 



AB 



^ AE^^AE'\bi^~ ACb, 



n^a 
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+didj 



-A-C+—E + 



Sail 



E - a^E" - 3aa' + 



/ 2 ^ 

na 



E' 



n 



/ aa a 
+2 — + 2- 



+ 



-A' - 2C' + 



naa 
\ a 



aa 



2d 



12 



n'aa' 



n 



1 • / n 
E--B+ - 



a 



n 
n 



A 



2n 



1 • a' ■ 



n 



an 

3a 



an^ 



{'271^ 2n^a 



B' 



an^ 



+ 2 



aa 



na 
n^a . 



B 



(1^ 77' (1^ 

- 2-A' - A' - 3C" - AS" - 2-AS', 

an a 



- + "-A 
n a 



a 



a 



+ \a-'AB' - ^AB 
2 a'^n 



, • , In' a'\ ^ ■ a , 

-AE' + - -A^;'. 

V n a / a 



B 



(45) 



(46) 
(47) 

(48) 



2. Vector components 



5Rl 



5Rl 



5Rl 



1 • 1 1 n' a'\ -, 
-AEi - -B' + R' + 

2 2* 2\n a ' 



a na 



a^^ 



+ 2- 



a'n' 



n^ 



1 1 ^' 

-ae; -S, 

2 * 2n2 * 



?7, 



3d 



3ad 



,2n'^ 2n^a 



a na 
B' + 



aa 



an 



an^ 



n^a^ 



Bi- 



na 
n^a 



duEi) - a^d^iE'^^ - 3aa' + 



/ 2 ^ 
na 



n 



d^.E'^^ 



+2 



aa 



a 



+ 2— - 



naa 



- aa" - 2a'2 - 



1 



ri a 

3 



n'aa' 



n 



d^iEj) 



an^ 



d(iBj) 



1 

2^ 
5, 



AR; 



(49) 
(50) 



(51) 



5. Tensor components 

ml = -1a£, + ^4 + ( ^ - ) 4 - 

ad a? naa 



/ 3ad 


a^n 




1 


1^2n2 


2n^ ^ 


2 


// 

— aa 


-2a'2 


n'aa' 


) 


n 






V. JUNCTION CONDITIONS 

In the previous section, we have considered only the perturbations in the bulk. It is now 
time to take into account the brane itself, and in particular the perturbations of the matter 
in the brane. The connection between the metric perturbations, living in the bulk, and the 



matter perturbations confined to the brane is made via the junction conditions |17 
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The energy-momentum tensor describing matter content of the brane will be written in 
the form 



= ^{y)S'., (53) 

and, in the background, it has the perfect fluid form 

S^^ = {p + P)u^u, + Pg^^. (54) 
The junction conditions are found to be given, in the five-dimensional case [Q, by 

[Kf,u\ = -K^ (s^^ - \sg^^ , (55) 



3 

where the brackets denote the jump across the brane, of the extrinsic curvature K^j_i, and 
S = Sj^. In a Gaussian normal coordinate system, the extrinsic curvature is given by the 
simple expression 

^M'^ = ^dyg^u- (56) 



In the present context, it is convenient to decompose the junction conditions ( ^5]) into 
a background part and a perturbed part. The use of the background junction conditions 
is necessary to find solutions of the Einstein equations when homogeneity and isotropy in 
the brane are assumed, such as those given in ( ^2] - ^ ). The junction conditions for the 
background have been given in M 



\a'] 



-^P, (57) 



aobo 3 

In'] , ^ , ^ 

^ = -3p + 2p, 58 
nobo 3 

where the subscript for a, h, n means that these functions are taken in ?/ = 0. 

Let us now consider the junction conditions for the perturbations, which can be written 

[SK^,] = ^ (^~6S^, + ^g^JS + i^v) , (59) 

Using the expression (^) for the extrinsic curvature in the GN gauge as well as the y —y 
symmetry, one ends up with the following condition, 

h'f,^\y=o+ = (-SS^u + ^g^^u^S + ^Sh^J^ . (60) 

It is then useful to decompose further these junction conditions by distinguishing space 
and time. Let us begin with the perturbations of the energy-momentum tensor in the brane. 
The perturbations for the unit four-velocity can be written 

Su^" = {-n-^A,a-\'} , (61) 
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where the expression for Su'^ follows automatically from the normalization condition satisfied 
by u'^. The perturbations of the energy-momentum tensor then have the following form: 

6Soo = n^Sp + 2pn^A, 

SSoi = -(p + P)navi - pBi, 

SSij = a^5P5ij + Phij + a^TVij, (62) 

with Vi = 5ijV^ and where TCij is the anisotropic stress tensor. Once more it is possible to 
decompose the above expressions for the brane matter energy-momentum tensor into scalar, 
vector and tensor components, using 

Vi = diV + Vi, (63) 

with div'' = and 



■Kij = l^didj - -^^ij^j + 29(i7rj) + TTy, (64) 

with as usual vTj transverse and vfjj transverse traceless. Substituting the above expressions 
for the perturbed energy-momentum tensor and using the background junction conditions 
( ^7| - |58|) yields the following conditions for the perturbations 

A\y=o^ = ^ {25 p + , (65) 
{n-'B;)\y^,^ = K\p + P)^v,, (66) 

h'ij\y=o+ = -1^^ (g-^P^ + T^ii) ■ (67) 

These conditions can also be decomposed, in a straightforward manner, into scalar, vector 
and tensor junction conditions. 



VI. CONCLUSION 

In the present work, we have developed a formalism in order to deal with the evolution 
of cosmological perturbations in a brane universe. This formalism has the advantage to 
introduce quantities that are similar to the usual treatments of cosmological perturbations 
in standard cosmology, thus making easier in the future a comparison between brane cos- 
moloical perturbations and observable quantities, such as the large scale structure and the 
temperature anisotropies. 

However, the equations governing the perturbations in the brane scenario are much 
more complicated than in standard cosmology. The main reason is that the equations for 
the metric perturbations, after a usual decomposition into ordinary (spatial) Fourier modes, 
will contain partial derivatives with respect to time and with respect to the fifth coordinate, 
in contrast with standard cosmology where one ends up with only ordinary differential 
equations with respect to time. Solving these equations appears to be quite a challenge. 

Note: while the present work was being completed, two works on the same subject ( |jT^ 
and |jl9|) have appeared on hep-th. 
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